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ABSTRACT 

The formation and evolution of the circumstellar disk in unmagnetized molec- 
ular clouds is investigated using three-dimensional hydrodynamic simulations 
from the prestellar core until the end of the main accretion phase. In collapsing 
cloud cores, the first (adiabatic) core with a size of ~ 10 AU forms prior to the 
formation of the protostar. At its formation, the first core has a thick disk-like 
structure, and is mainly supported by the thermal pressure. After the protostar 
formation, it decreases the thickness gradually, and becomes supported by the 
centrifugal force. We found that the first core is a precursor of the circumstellar 
disk. This indicates that the circumstellar disk is formed before the protostar 
formation with a size of ~ 10 AU, which means that no protoplanetary disk 
smaller than < 10 AU exists. Reflecting the thermodynamics of the collapsing 
gas, at the protostar formation epoch, the first core (or the circumstellar disk) 
has a mass of ~ 0.01 — 0.1 M & , while the protostar has a mass of ~ 1O _3 M . 
Thus, just after the protostar formation, the circumstellar disk is about 10 — 100 
times more massive than the protostar. Even in the main accretion phase that 
lasts for ~ 10 5 yr, the circumstellar disk mass dominates the protostellar mass. 
Such a massive disk is unstable to gravitational instability, and tends to show 
fragmentation. Our calculations indicate that the planet or brown-dwarf mass 
object may form in the circumstellar disk in the main accretion phase. In ad- 
dition, the mass accretion rate onto the protostar shows strong time variability 
that is caused by the perturbation of proto-planets and/or the spiral arms in the 
circumstellar disk. Such variability provides a useful signature for detecting the 
planet-sized companion in the circumstellar disk around very young protostars. 
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Introduction 



We believe that stars are born with a circumstellar disk. The formation of the circum- 
stellar disk is coupled with "the angular momentum problem" that is a serious problem in 
the star formation process, and the dynamics of disks may determine the mass accretion rate 
onto the protostar that determines the final stellar mass. In addition, planets are considered 
to form in the circumstellar (or protoplanetary) disk, and their formation process strongly 
depends on disk properties such as disk size and mass. Thus, the formation and evolution 
of the circumstellar disk can provide a significant clue to star and planet formation. 



Stars form in molecular clouds that h ave an angular momentum (jArquilla fc Goldsmith 



19861 : iGoodman et al. Ill993l ; ICaselrill2002l ). Thus, the appearance of a circumstellar disk is a 
natural consequence of the star formation process when the angular momentum is conserved 
in the collapsing cloud core. In additio n, observations hav e shown the existence of circum- 
stella r disks around the protostar (e.g., IWatson et al.l 120071 ; iDutrey et al.l 120071 : iMeyer et al. 



20071 ). Numerous observations indicate that circumstellar disks around C lass I and II proto 



stars have a size o f ~ 10- 1000 AU and a mass of ~ 1O" 3 -O.1M (e.g., ICalvet et al.ll2000 



Natta et al.l 120001 ) . However, they correspond to phases long after their formation. Because 
the formation site of the circumstellar disk and protostar are embedded in a dense infalling 
envelope, we cannot directly observe newborn or very young circumstellar disks (and proto- 
stars). Thus, we only observe the circumstellar disks long after their formation, i.e., around 
the class I or II protostar phase. 

Observations also indicate that a younger protostar ha s a m assive circumstellar disk 
( INatta et al.ll2000l : IMeyer et al.ll2007l ). Recently, Enoch et al.l ( 120091 ) observed a massive disk 
with M(j is k ~ 1 M Q around class sources, indicating that this massive disk can be present 
early in the main accretion phase. However, unfortunately, observation cannot determine 
the real sizes of circumstellar disks, and how and when they form. Therefore, we cannot 
understand the formation process of the circumstellar disk by observations. The theoretical 
approach and numerical simulation are necessary to investigate the formation and evolution 
of the circumstellar disk. 

Theoretically, the star formation process can be divided into two phases, i.e., the early 
collapse phase and main accretion phase. The molecular cloud cores that are star cradles 
have a number density of n ~ 10 4 — 10 6 cm~ 3 . In the early collapse phase, the gas in the 



-3 - 



molecular cloud core continues to collapse until the protostar formation at n ~ 10 cm . 
The p rotostar at its formation has a mass of M p s ~ 10~ 3 M that corresponds to the Jovian 



mass (lLarsonlll969l ; lMasunaga &: Inutsuka Il2000l ). In the main accretion phase, the protostar 
acquires almost all its mass by the gas accretion to reach ~ 1M . In this paper, we define 
'the early collapse phase' as the period before the protostar formation, while 'the main 
accretion phase' is defined as the period until the gas accretion onto the protostellar system 
(protostar and circumstellar disk) almost halts after the protostar formation. In general, 
it is considered that the circumstellar disk gradually increases its mass and size in the 
main accretion phase that is successively connected from the early collapse phase. Thus, to 
understand the formation and early evolution of circumstellar disks, we should consider both 
the early collapse and main accretion phases; we need a self-consistent calculation from the 
collapse of the molecular cloud core until the end of the main accretion phase through the 
protostar formation. 

However, to investigate the formation and evolution of circumstellar disks in numerical 
simulations, we need a very long-term calculation with a sufficient spatial resolution, in which 
we should calculate the evolution of the protostellar system at least for the time comparable 
to the freefall timescale of the initial cloud core, i.e., ^ 10 4 — 10 5 yr after the protostar 
formation. In addition, we should resolve spatial scale length down to at least ~ 1AU. 
Reflecting the thermodynamics of the coll apsing gas, two nested cores with a typical spatial 



scale appear in the early collapse phase ( IL arson! Il969l ; iMasunaga fc Inutsuka 



20001 ). The 



inner core (so-called the second adiabatic core) corresponds to the protostar that has a 
size of ~ 1-R©, while the outer cor e that is called the fir st (adiabatic) core has a size of 
~ 1 — 10 AU (for details, see §3.1. ip . Ilnutsuka et al. I (120091 ) expected that the circumstellar 
disk originates in the first core formed in the early collapse phase. Thus, to investigate the 
formation of the circumstellar disk, we have to resolve the first core spatially; we need a 
spatial resolution of <^ 1 AU at least. Moreover, we need a three dimensional calculation to 
properly treat the angular momentum transport in the circumstellar disk. 

So far, many studies o f star formation in the collapsing cloud mainly foc used only on the 



early collapse phase (e.g., iBodenheimer et al. Il2000l ; iGoodwin et al. 1 120071 ). To investigate 



the formation and evolution of the circumstellar disk, we have to calculate the early col- 
lapse and subsequent main accretion phases. However, since such calculation requires a huge 
amount of CPU time, only a few studies reported the formatio n of the circumstella r disk in 
the collapsing cloud core including the main accretion phase. Kratter et al.l (120091 ) investi- 
gated the formation of the circumstellar disk with sub-AU resolut ion under the isoth ermal 
approximation, and showed frequent fragmentation of the disk. IWalch et al.l (120091 ) also 
studied the circumstellar disk formation in the collapsing cloud core with a slightly coarser 
spatial resolution of 2AU approximating the radiative cooling with adiabatic equation of 
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state, and showed properties of the circumstellar disk. 



In this study, using three-dimensional simulation adopting a barotropic equation of 
state with a higher spatial resolution than in previous studies, the evolution of the unmag- 
netized collapsing cloud cores is investigated until the gas accretion onto the protostar and 
circumstellar disk almost halts. In three dimensions, we first calculate the evolution of the 
circumstellar disk by the end of the main accretion phase. Although the circumstellar disk 
formation in the collapsing cloud may be investigated with a radiation-hydrodynamics code, 
it is very difficult to execute such a calculation even with current supercomputers because 
it takes a huge amount of CPU time. Thus, to trace the gas thermodynamics, we chose a 
barotropic equation of state that is used even in recent two-dimensional simulations of the 
circumstellar disk formation (e.g., Vorobyov fc Basu 2007 ; Vorobyovl 2009 ). The structure of 



the paper is as follows. The framework of our models and the numerical method are given 
in §2. The numerical results are presented in §3. We discuss the fragmentation condition 
of the circumstellar disk and implication of the planet formation in §4, and summarize our 
results in 85. 



2. Model Settings 



To study the evolution of collapsing gas clouds and circumstellar disks, we solve the 
equations of hydrodynamics including self-gravity: 



dp 
dt 



+ v ■ (H = o, 



V 2 6 = 



= -VP 

4-KGp, 



pV0, 



(2) 
(3) 



where p, v, P, and <fi denote the density, velocity, pressure, and gravitational poten t ial, re - 
spectively. To mimic the temperature evolution calculated by iMasunaga fc Inutsuka I ( 120001 ). 
we adopt the piece-wise polytropic equation of state (see, iMachida et al.ll2007h as 



c loP 



where c Sj o = 190 ms 1 , and p c = 3.84 x 10 14 g cm 
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As the initial state, we take a spherical cloud with critical Bonnor-Ebert (BE) density 
profile, in which the uniform density is adopted outside the sphere (r > R c ). For the BE den- 
sity profile, we adopt the central density of n c = 6 x 10 5 cm -3 and isothermal temperature of 
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T = 10 K. For these parameters, the critical BE radius is R c = 6.2 x 10 3 AU. To promote the 
contraction, we increase the density by a factor of /=1.68, where / is the density enhance- 
ment factor that represents the stability of the initial cloud. With / = 1.68, the initial cloud 
has (negative) gravitational energy twice that of thermal energy. Thus, the central density 
of the initial sphere is n C) ; n i = 10 6 cm -3 , while the ambient density is n am b = 7.2 x 10 4 cm~ 3 . 
We add m = 2-mode non-axisymmetric density perturbation to the initial core. Then, the 
density profile of the core is described as 

n(r ) = / Pbe (x) (1 + S„) f for r < R c , 

P{ 1 I P BE(Rc)(l + 5 p )f for r>R Cl 1 1 

where Pbe(^) is the density distribution of the critical BE sphere, and S p is the axisymmetric 
density perturbation. For the m = 2-mode, we chose 

5 P = Mr/Rc) 2 cos 20, (6) 

where (=0.01) represents the amplitude of the perturbation. The radial dependence is 
chosen so that the density perturbation remains regular at the origin (r = 0) at one time- 
step after the initial stage. This perturbation ensures that the center of the gravity is always 
located at the origin. The mass within r < R c is M = 1M . The gravitational force is 
ignored outside the host cloud (r > R c ) to mimic a stationary interstellar medium. Initially, 
the cloud rotates rigidly with angular velocity Qq around the z-axis. We parameterized 
the ratio of the rotational to the gravitational energy (/3 ) inside the initial cloud. With 
different /3 , we calculated 10 models. Model names, initial angular velocities fl , and (3 are 
summarized in Table 1. 

In the collapsing cloud core, we assume protostar formation occurs when the number 
density exceeds n > 10 13 cm -3 at the cloud center. To model the protostar, we adopt a sink 
around the center of the computational domain. In the region r < r sin k = 1 AU, gas having 
a number density of n > 10 13 cm -3 is removed from the comp utational domain and added 



to the protostar as a gravity in each timestep (for details, see iMachida et al.ll2009l ). This 
treatment of the sink makes it possible to calculate the evolution of the collapsing cloud and 
circumstellar disk for a longer duration. 

To calculate oyer a large spa tial scale, the nested grid method is adopted (for details, 



sec 



Machida et al. Il2005al . l2006al ). Each level of a rectangular grid has the same number of 
cells of 128 x 128 x 16. The calculation is first performed with five grid levels (I = 1-5). The 
box size of the coarsest grid I = 1 is chosen to be 2 5 R C . Thus, a grid of I = 1 has a box size 
of 1.97 x 10 5 AU. A new finer grid is generated before the Jeans condition is violated. The 
maximum level of grids is restricted to Z max = 12. The / = 12 grid has a box size of 96 AU 
and cell width of 0.75 AU. With this method, we cover five orders of magnitude in spatial 
scale. 
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Results 



We investigated the cloud evolution and the circumstellar disk formation with different 
initial rotational energies, /3q. Observations have sh own that the molecular cloud cores have 
10" 4 £ 0o £ 0.02 with the mean value of A, ~ 0.02 jGoodman et al. Ill993l : ICasellill2002h . In 
the following section, we show the cloud evolution and formation process of the circumstellar 
disk for two typical models (model 6 and 3). Then, we compare the properties of circumstellar 
disks in clouds with different rotational energies. 



3.1. Typical Models 

3.1.1. Cloud with 0o = 10~ 3 

Figures [1] and [2] show the time sequence around the center of the cloud before (panels a 
and b) and after (panels c - f) the protostar formation for the model with O = 10 -3 (model 
6). As denoted in §2J we define the protostar formation epoch (t c = 0) as the time when 
the maximum density reaches n = 10 13 cm -3 . In these figures, the time after the protostar 
formation (t c ) and protostellar mass (M ps ) are described in each panel. The protostellar 
mass M ps is derived as the mass falling into the sink. 

Panels a and b in these figures indicate the formation of the disk-like structure before 
the protostar formation. To clearly define the disk, we used the ratio of the radial to 
azimuthal velocity R v (= |i> r /ity|). Figure [3] shows the density distribution (a), the ratio of 
the radial to azimuthal velocity (6), and the distribution of the radial (c) and azimuthal (d) 
velocities against the distance from the protostar along the y-axis. In Figure [31 these values 
for four different epochs corresponding to panels b, c, e and / of Figures [T] and [2] are plotted. 
Comparison of Figures [1] and [2] with Figure [3] a indicates that the shock front corresponds 
to the disk surface. Thus, the disk-like structure is enclosed by the shock. In Figure [3]a, fine 
structures inside the shock surface at t c = 2.1 X 10 4 yr and at t c = 1.0 x 10 5 yr correspond to 
spiral density waves. Figure [3] (c) shows a sudden rise of the radial velocity at the shock front. 
In addition, Figure [3]d shows a large azimuthal velocity inside or near the shock front (i.e., 
inside or near the disk). Figures [3]c and d indicate that the gas rapidly falls into the center of 
the cloud with a slow rotation outside the disk, while the radial velocity slows and azimuthal 
velocity dominates inside the disk. Therefore, the ratio of the radial to azimuthal velocity 
shows a sudden drop at the disk surface as shown in Figure Note that, in Figured, to 
stress the disk surface, the velocity ratio R v inside the disk is not displayed. As a result, the 
velocity ratio R v is a good indicator to specify the disk. 
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In this paper, to determine the disk, we estimated the velocity ratio R v (= v r jv^) in 
each cell, and specified the most distant cell having R v < 1 from the center of the cloud. 
Then, we defined the disk radius r disk as the distance of the cell furthest from the origin, and 
the disk boundary density pd,b as the density of the most distant cell. Finally, we defined 
the disk that has a density of p > pd,b- Thus, in our definition, the region inside the rapid 
drop of R v corresponds to the disk in Figure [3]&. Comparison of Figure |3f> with Figures [fl 
|2] and [3]a indicates this definition of the disk well corresponds to a real size of the disk in 
simulation. 

Figure [2] shows that a thick disk-like structure formed before the protostar formation 
(Fig. |2] a and b) transforms into a sufficiently thin disk after the protostar formation (Fig. [2] 
c — f). As described in theoretically, the star formation process can be divided into two 
phases: the early collapse phase (or the early phase of the star formation) and the main 
accretion phase (or the later phase of the star formation). The early c ollapse phase was 



investigated in detail by iLarsonl ( 119691 ) and iMasunaga fc Inutsuka I (120001 ) . Here, we briefly 



describe this. After the gas collapse is initiated in the molecular cloud core, the collapsing 
gas obeys the isothermal equation of state with temperature of ~ 10K for n c ^ 10 10 cm -3 
(isothermal phase). Then cloud collapses adiabatically (10 10 £ n c <^ 10 16 ; adiabatic phase) 
and quasi-static core (i.e., first core) forms during the adiabatic phase. After central density 
reaches n c ~ 10 16 cm -3 , the equation of state becomes soft reflecting the dissociation of 
hydrogen molecules at T ~ 2 x 10 3 K, and the gas collapses rapidly again, i.e., the second 
collapse begins. Finally, when the gas density reaches n c ~ 10 21 cm -3 , the gas collapse stops 
and the protostar (or the second core) forms. At this epoch, the ea rly collapse phase ends 



and th e main accretion phase begins. In the main accretion phase, IMasunaga fc Inutsuka 



(120001 ) expected that the first core without the angular mom entum disappears in ~ 100 yr 



after the protostar formation, while ISaigo fc Tomisakal (120061 ) pointed out that the first core 



having the angular momentum does not disappear in such short duration because the first 
core is supported by the rotation. Thus, we can expect that the rotating first core formed 
in the early collapse phase becomes the circumstellar disk in the main accretion phase. In 
the main accretion phase, the gas with the angular momentum continues to accrete onto the 
first core (or the circumstellar disk). 

Figure [2] clearly shows that the first core becomes the circumstellar disk in the main 
accretion phase; the first core is a precursor of the circumstellar disk. The region enclosed 
by the shock in Figure [T]a and Figure [2]a corresponds to the first core that is formed after 
the gas becomes adiabatic. Figures [2]&- e show that the first core (or the circumstellar disk) 
gradually becomes thin with time. Finally, a sufficiently thin disk appears as shown in 
Figure [2jf. Figure [3]& shows that the disk has a radius of ~ 10 AU before the protostar 
formation. The disk extends up to ~ 200 AU for t c ~ 10 5 yr in the main accretion phase. In 
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addition, 89% of the total mass (i.e., 89% of the initial host cloud mass) accretes onto the 
protostar and circumstellar disk system by this epoch (t c <^ 10 5 yr, see Table 1). 

In summary, when we defined the circumstellar disk as the rotating disk around the 
protostar, the circumstellar disk with a size of ~ 10 — 20 AU is already formed before the 
protostar formation. T his disk size at the protostar formation corresponds to that of the 
first core with rotation (jMatsumoto fc Hanawal 120031 : ISaigo &: Tomisakall2006l ). Therefore, 
the circumstellar disk has a minimum size of ~ 10 AU which implies that we cannot observe 
a disk with ^ 10 AU even around very young protostars. 

As shown in Figure (U fragmentation occurs and two clumps form in the circumstellar 
disk ~ 1.5 x 10 4 yr after the protostar formation. At the fragmentation epoch, the protostellar 
mass is M ps ~ 0.1 M , while the disk mass is M^w ~ 0.18 M . Thus, the disk is 1.8 ti mes 



more massive than the protostar. As shown in iLarsonl (jl969l ) and lMasunaga et al. I (jl998l ). at 
the protostar formation epoch, the protostar only has a mass of ~ 10 -3 M & that corresponds 
to Jovian mass. This Jovian mass protostar acquires its mass by the gas accretion for 
~ 10 5 — 10 6 yr to re ach the solar mass object. On the other hand, the first core has a mass 
of ~ 0.01 — 0.1 M & (jMatsumoto fc Hanawal 120031 ). The mass of the protostar and first core 



correspond to the Jeans mass at their formation epoch. Figures [TH3] show that the first core 
evolves directly into the circumstellar disk through the early collapse to the main accretion 
phases. Thus, during the early main accretion phase, the ci rcumstellar disk is more massive 



than the protostar, as described in llnutsuka et al. I (120091). Since such a massive disk is 



gravitationally unstable, fragmentation tends to occur (jToomrelll964j ). In this model, at the 
end of the calculation, each fragment is gravitationally bound, and has a mass ~ 0.037 M & 
with ~ 12 AU of the separation between the central protostar and fragment. 



3.1.2. Cloud with p = 10~ 2 

Figure H] shows the time sequence for the model with /3 = 10 -2 (model 3). Similar 
to the model with (3 = 10~ 3 (model 6), this model also shows the disk-like structure (i.e., 
disk-like first core) preceding the protostar formation (Fig. H]a). Owing to the larger initial 
angular momentum, the first core in model 3 (/3o = 10~ 2 ) has a flatter structure than that 
in model 6 (/?o = 10~ 3 ). For this model, at the protostar formation epoch, the protostar 
has a mass of M ps ~ 10~ 3 M Q , while the disk has a mass of Mdi S k — 0.1 M Q . Thus, the disk 
is about 100 times more massive than the protostar. Nevertheless, the circumstellar disk 
shows no fragmentation, as seen in Figure H]& - d. Instead, the spiral structures appear with 
the gravitational instability in the circumstellar disk. The circumstellar disk extends up to 
~ 500 AU (Fig. 0i) by the epoch at which the almost all gas in the host cloud has accreted 
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onto the protostar and circumstellar disk system. By this epoch, the masses of protostar and 
circumstellar disk reach M ps = 0.14 M and M disk = 0.58 M , respectively. Thus, even at 
end of the main accretion phase, the circumstellar disk is more massive than the protostar. 
The disk properties and fragmentation condition are discussed in §4.11 

3.2. Disk Properties vs. Initial Rotational Energies 

In this section, to investigate the disk properties, four different models with different 
initial rotational energies (3o are presented. Figure [5] shows the density and velocity distri- 
bution for models with /3q = 10 -2 , 10~ 3 , 10~ 4 and 10~ 5 at t c ~ 10 5 yr after the protostar 
formation. Since the freefall timescale at the center of initial cloud is tff — 2.3 x 10 4 yr, 
they are structures at ~ 4.3 tg t a after the protostar formation. The figure shows a larger disk 
with larger (3q. In addition, two clumps formed by fragmentation due to the gravitational 
instability appear in the circumstellar disk for models with (3 = 10~ 3 and 10~ 4 , while no 
clump appears by the end of the main accretion phase for models with O = 10 -2 and 10~ 5 . 
The fragmentation condition depends on the size and mass of the circumstellar disk (see, 

Figure M shows the evolution of the disk mass before (left panel) and after (right panel) 
the protostar formation for the same models in Figure The figure indicates that the 
rotating disk forms before the protostar formation (i.e., t c < 0) and has a mass of 6 x 10~ 3 — 
0.1 M & at the protostar formation epoch. After the protostar forms, the circumstellar disk 
gradually increases its mass by gas accretion, and reaches M disk = 0.1 — 0.6 M by the end 
of the main accretion phase. The mass accretion rate onto the circumstellar disk during the 
main accretion phase (~ 10 5 yr) is M d i S k = (1 — 5) x 10 -6 M yr _1 . This rate is determined 
by the balance between the mass accreting onto the circumstellar disk and mass infalling 
onto the protostar. 

Figure [7]a plots a time sequence of the residual mass M res = M ini — M disk — M ps for 
the same models in Figure where M ini is the initial mass of the host cloud (i.e., the mass 
inside r < R c in the initial cloud). The figure shows that the residual mass rapidly decrease 
in a freefall timescale and reaches ~ 0.1 M Q in ~ 10 5 yr. Thus, about 90% of total mass 
accretes onto the protostar and circumstellar disk system by this epoch (see also, Table 1). 
Therefore, the gas accretion almost halts and the main accretion phase ends at t c ~ 10 5 yr. 

Figure [7]6 shows the time evolution of the protostellar mass for the same models. The 
protostar has a mass of (1 — 3) x 10 -3 M Q at its formation epoch. Then, the protostar 
acquires its mass by gas accretion and reaches 0.3 — 0.9 M by the end of the main accretion 
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phase. The mass accretion rate onto the protostar during the main accretion phase is M ps = 
(3 — 9) x 10~ 6 M Q yr -1 on average. The time variability of the mass accretion rate is shown 
in §3.31 Figures O and [7]& indicate that the model with a larger rotational energy has a larger 
Miisk, but smaller M ps . As shown in Figure 0, the cloud with a larger (3 has a larger (or 
massive) rotating disk. In such a disk, the mass accretion onto the protostar is suppressed 
owing to larger angular momentum, and a less massive protostar appears with a relatively 
smaller protostellar accretion rate. 

Figure M shows the disk-to-protostellar mass ratio \i against the protostellar mass for 
the same models in Figure [5j The figure indicates that the circumstellar disk mass exceeds 
protostellar mass at the protostar formation epoch even for the model with a considerable 
small rotational energy. The cloud with larger rotational energy has a relatively massive 
circumstellar disk at the protostar formation epoch. For example, the circumstellar disk for 
model with (3q = 10 -2 is about 100 times more massive than the protostar (/i = 100), while 
for the model with /3 = 10~ 5 it is twice as massive (/i = 2). The circumstellar disk is more 
massive than the protostar for most of the main accretion phase, exception for the model 
with /?o = 10~ 5 . In addition, the circumstellar disk mass is comparable to the protostellar 
mass even at the end of the main accretion phase for models with /3 > 10~ 4 . Moreover, even 
in a cloud with /3n = 10 ~ 5 whose value is much smaller than the average of the observation 
0o = 0.02 (jCasellill2002l ). the circumstellar disk has a 1/10 (/i = 0.1) of the protostellar mass 
at the end of the main accretion phase. 

Figure He plots the evolution of the circumstellar disk radius after the protostar forma- 
tion. The figure shows that the circumstellar disk with the radius of 3 — 40 AU already exists 
at the protostar formation epoch (i.e., t c = 0). The size of the circumstellar disk increases 
with time, and reaches 30 — 500 AU by the end of the main accretion phase. The figure also 
indicates that the protostar has a larger circumstellar disk for the cloud with larger j3o- 

Figure Wd plots the aspect ratio of the circumstellar disk, H/R, that is defined as 
H/R = h z /(hih s Y^ 2 , where hi, h s , and h z are, respectively, lengths of the m ajor, minor, and 



z-axes derived from the moment of inertia inside the circumstellar disk (see, Machida et al. 



2005aj). Note that since hi, h s and h z are weighted by the mass, they are slightly different 



from a real size of the circumstellar disk which is surrounded by the accretion shock as shown 
in Figures [1] and |2j The figure shows that, for models with /?o < 10~ 2 , the protostars have 
relatively thick circumstellar disks with H/R > 0.5 at their formation epochs as also seen 
in Figure [2]c. On the other hand, for model with /3 = 10~ 2 , the protostar at its formation 
already has a thin circumstellar disk with H/R < 0.2. However, in any model, the aspect 
ratio decreases with time and reaches H/R < 0.1 by the end of the main accretion phase. 
Thus, protostars have a sufficiently thin disk at the end of the main accretion phase. 
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3.3. Mass Accretion Rate onto Protostar 

Figure |9] shows the mass accretion rate onto the protostar M ps (left axis) and circum- 
stellar disk mass (right axis) against the protostellar mass. As described in ^21 we removed 
the gas having a number density of n > 10 13 cm -3 in the region r < 1 AU from the com- 
putational domain. We regard the removed gas as the accreting mass onto the protostar, 
and estimated the mass accretion rate in each step. The figure shows that the accretion 
rates are in the range of 10~ 4 ^ M ps /(M yr -1 ) <^ 10 -6 in each model. In the theoretical 
analysis of star formation, the mass accretion r ate is prescri bed as M = f c 3 /G, where / is a 



numerical factor (e.g., / = 0.975 for IShulll977l . / = 46.9 for iHunterl 119771 ) . Since gas clouds 
have temperatures of T = 10 K (c s = 0.2 km s" 1 ), the accretion rate in the main accretion 
phase is M = (2 — 90) x lO~ 6 M yr . Thus, the accretion rate derived in our calculations 
corresponds well to the theoretical expectation. However, the accretion rate shows strong 
time variability as the protostellar and circumstellar mass increases. 

As seen in Figure EE for models with {3 < 10~ 2 , the accretion rate monotonically 
decreases initially, while it strongly fluctuates with time later as the protostellar (or circum- 
stellar) mass increases. In addition, models with /3q = 10~ 2 show a strong time variability 
of the accretion rate right after the protostar formation. This time variability stems from 
spiral arms or clumps caused by the gravitational instability in the circumstellar disk. The 
epoch at which the accretion rate begins to fluctuate corresponds to the formation epoch of 
the spiral arm or clump. For example, for model with /3q = 10~ 3 , the accretion rate shows a 
slight time variability when M ps ^ 0.04 M & (see, Fig. [T]c and d), while it shows a strong vari- 
ability when M ps ^ 0.04 M (see, Fig{T]e-/). Indeed, when M ps ~ 0.04 M , snapshot shows 
spiral structure. Thus, this time variability of the mass accretion rate onto the protostar is 
closely related to the stability of the disk. We discuss the disk stability in §4.11 



4. Discussion 

4.1. Stability and Fragmentation Condition for Circumstellar Disks 

As shown in £J3], the circumstellar disk has a mass larger than or comparable to the 
protostellar mass in the main accretion phase. Such massive disk tends to show fragmentation 
or the formation of spiral structure that is closely related to the mass accretion onto the 
central star, subsequent evolution of the circumstellar and protoplanetary disk, and planet 
formation. In this subsection, we discuss the disk stability and fragmentation condition. 

As seen in Figure [5j for typical models (models 3, 6, 8, and 10), fragmentation occurs 
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in clouds with moderate rotational energy of /3q = 10~ 3 and 10~ 4 , while spiral arms appear 
without fragmentation in model with larger (/3 = 10 -2 ) and smallest (/3q = 10 -5 ) rotational 
energies. T o investigate t he fragmentation condition of the disk, we calculated Toomre's Q 
parameter ( IToomrel 119641 ). which is defined as 



Q 



C S K 



(7) 



where c s , k and £ are the sound speed, the epicyclic frequency, and the surface density, 
res pectively. We also calculated the ratio of the critical Jeans length A cr j to the disk radius 
r<j ( IMatsumoto fc Hanawal 120031 1 as 

R c 



Acri 



where the critical Jeans length can be described as 



A, 



[i + (i - Q) 1 ' 2 ] 



(9) 



It is considered that fragmentation occurs when both conditions of Q < 1 and R c < 1 are 
fulfilled in the disk. At eac h timestep, we estimated Q a nd A c in each cell inside the disk 
(for detailed description, see IMatsumoto fc Hanawal 120031 ) . Then, these values are averaged 
within the whole disk to roughly estimate the time variability of them. Note that although 
we may have to check local values of Q and R r , not global (or averaged) values to investigate 
the disk stability and fragmentation (IKratter et al.ll2009l ). we use global values to roughly 
estimate the time sequence of the disk stability. 

As shown in £j3j the circumstellar disk originates from the first core that is formed in 
the early collapse phase. The first core is formed after the gas becomes adiabatic in the 
collapsing cloud core. Thus, the first core is mainly supported by the thermal pressure 
gradient force at its formation. Even when the initial host cloud has a significant rotational 
energy, the rotational energy of the first core is, at most, only 10% of the gravitational energy 



(IMachida et al. Il2005al l2006al ). Then, the gas having a larger angular momentum accretes 



onto the disk (the first core or circumstellar disk) from the outer envelope as time goes on, 
and the disk gradually thins owing to the centrifugal force, and spreads to a large extent. 
Thus, it is considered that Q increases with time, because the gas with the larger specific 
angular momentum accretes onto the circumstellar disk with time; the epicyclic frequency 
or the ratio k/H is expected to become large with time. Figure [10] left panel plots a time 
sequence of the averaged Q. The figure shows that, for any model, Q is much less than 
unity (Q <C 1) at the protostar formation epoch and increases with time (or the protostellar 
mass). Finally, it reaches order unity Q ~ 1 by the end of the main accretion phase in 
any model. The figure also shows that the circumstellar disk has a smaller Q in the cloud 
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with smaller rotational energy. This is because, for the model with smaller ft, the gas with 
smaller specific angular momentum forms a disk. 

In general, fragmentation can occur in the rotating disk with Q < 1. However, for the 
model with ft = 10 -2 , although the circumstellar disk has Q < 1 when the protostellar 
mass is less than M ps < 0.02 M , fragmentation does not occur, as shown in Figure HI 
This is because the disk radius is smaller than the critical Jeans length (R c > 1) when 
M ps < 0.02 Mq. On the other hand, for this model, the disk radius is larger than the critical 
Jeans length R c < 1 when protostellar mass is larger than M ps > 0.4, while fragmentation 
does not occur because Q parameter exceeds unity Q > 1 in this phase. As a result, by the 
end of main accretion phase, this model shows no fragmentation, while the spiral structure 
appears as shown in Figure HI 

For the model with ft = 10 -3 , the circumstellar disk keeps Q < 1 for M ^ 0.2 M , and 
the ratio R c manages to reach R c < 1 in this period. Therefore, fragmentation occurs to 
form two clumps as shown in Figures [I] and [2j Even in the model with ft = 10 -4 , these two 
conditions (Q < 1 and R c < 1) are fulfilled for M ps <^ 0.1 M , and fragmentation occurs. 
On the other hand, the model with ft = 10~ 5 has Q < 1 through the main accretion phase, 
while fragmentation does not occur. This is because the disk size is smaller than the critical 
Jeans length by the end of the main accretion phase. 

To investigate the fragmentation condition in the main accretion phase, we calculated 
10 models in total, as listed in Table 1. Figure O implies that the upper border of the 
rotational energy between fragmentation and non-fragmentation models exists in models 
between ft = 10~ 2 and 10~ 3 . To determine the fragmentation condition in more detail, we 
show the final fate for models ft = 3 x 10~ 3 , 5 x 10~ 3 and 3 x 10~ 2 in Figure [TTJ The 
figure indicates that fragmentation border exists in models between ft = 3 x 10~ 3 and 
5 x 10~ 3 . Fragmentation occurs in the model with ft = 3 x 10~ 3 , while fragmentation does 
not occur in the model with ft = 5 x 10~ 3 . In addition, fragmentation occurs again in the 
model with ft = 3 x 10~ 2 . Note that, for this model, we stopped calculation long before 
the end of main accretion phase (t c = 2.8 x 10 3 yr) because fragments move outwardly and 
the Jeans criterion is violated. This model shows fragmentation in the early main accretion 
phase of t c = 4 x 10 3 yr (see, Table 1), because a sufficiently large disk is already formed in 
the early collapse phase (i.e., before the main accretion phase). Moreover, the model with 
ft = 5 x 10~ 2 shows fragmentation in the early collapse phase (i.e., before the protostar 
formation, see, Table 1). T he fragmentation process in the early collapse pha se was reported 



in many past studies (e.g., iBodenheimer et al. Il2000i ; iGoodwin et al. 1 120071 ). In summary 



the fragmentation condition derived in our simulation can be described as follows: 
1. ft > 3 x 10~ 2 : fragmentation in the early collapse phase, 
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2. 10 5 < 0o < 5 x 10 3 : fragmentation in the main accretion phase, 

3. 0o < 10~ 5 and 5 x 10~ 3 < (3 < 3 x 10~ 2 : no fragmentation. 

Although models having parameters of 5 x 10~ 3 < 0o < 3 x 10~ 2 show no fragmentation, these 
models may show fragmentation in subsequent evolution stages. As listed in Table 1, the 
mass of the circumstellar disk for no fragmentation models is comparable to the protostellar 
mass at the end of the main accretion phase, in which Toomre's Q is nearly unity Q ~ 1. 
Such a disk can fragment to form several clumps when the cool ing time scale is shorter 
than the Kepler rotational timescale in the protoplanetary disk ( Gammiel boOll ). In this 



study, we adopted a barotropic equation of state that makes the direct disk calculation 
from the molecular cloud until the end of the main accretion phase possible. In the main 
accretion phase, the mass of the circumstellar disk overwhelms the protostellar mass with 
Q < 1. In such a situation, it is expected that the cooling process is not so important for 
fragmentation. Ins tead, to investigate a marginal stable disk with Q ~ 1, the cooling process 



may be important (IDurisen et al.ll2007l ). Ideally, we should investigate the circumstellar disk 
formation from the molecular cloud core with a radiative hydrodynamics, while we need to 
wait for further development of computational power to perform such calculations. However, 
we can expect that fragmentation in the circumstellar disk generally occurs because a massive 
circumstellar disk inevitably appears in the star formation process. 



4.2. Implication for Planet Formation 



As described in §4.11 a massive circumstellar disk inevitably appears in the main ac- 
cretion phase. This is because the first core that is a precursor of the circumstellar disk is 
about 10 — 100 times more massive than the protostar, reflecting the thermodynamics of the 
collapsing g as. Such a massive disk tends to show fragmentation owing to the gravitational 
instability (IDurisen et al.l 120071 ). In our calculation, many models showed fragmentation 
and formation of small clumps in the circumstellar disk. As listed in Table 1, at the end 
of the calculation, fragments have a mass of Mf rag = (19 — 72) x 10 _3 M Q with a separa- 
tion of r sep = 7.1 — 94.2 AU. For every model, fragments survived without falling into the 
central protostar. Although the masses of fragments exist in the range of the brown dwarf 
(Mf rag = 13 — 75 Mj up ), the planet-mass object may appear in the subsequent star formation 
phase, in which it is expected that the circumstellar disk cools to form less massive clumps 
(IDurisen et al.l 120071 ) . In this study, we calculated the evolution of the circumstellar disk 
until the main accretion phase ends from the prestellar core phase. However, further long- 
term calculations including an adequate radiative effect with a higher spatial resolution are 
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necessary to properly investigate the mass and separation of fragments. Since such calcula- 
tion is beyond the scope of this study, we qualitatively, not quantitatively, discuss the planet 
formation in the subsequent evolution phase of the circumstellar disk below. 



Recently, brown dwarf or planet mass companions have been observed. iNeuhauser et al. 



(120051 ) observed a brown-dwarf mass companion around GQ Lup with sepa ration of 



100 A U. In 2008, di rect images of plan et mass companions around Fomalhaut ( iKalas et al. 



20081 ) and HR8799 (IMarois et al.ll2008l ) are shown, in which they have masses of 3 — 10 Mj up 
with separations of 24—98 AU. It is considered that there are two different ways to form plan- 
ets in the circumstellar (or protoplanetary) disk: one is t he core accretion model ( IHayashi et al 
19851 ). the other is the gravitational instability model (jCameronl Il978l ) . In the core accre- 
tion scenario, it is difficult to account for the formation of massive planet in the region far 
from the protostar (r sep >, 5 — 10 AU) in the standard disk model ( iPerri fc Cameron! Il974 



Mizuno et al.l Il978t IHayashi et al. Ill985l ). On the other hand, a massive planet can form 
even in the region far from the protostar by fragmentation due to gravitational instability 
when the disk is sufficiently massive. In this study, we showed that, in the circumstellar disk, 
gravitationally bound clumps can form in the region of r sep > 5 AU. In addition, our calcula- 
tion indicates that the massive circumstellar disks comparable to or more massive than the 
protostar form in the main accretion phase when the r otational ene rgy of the molecular cloud 
is comparable to the observations (10~ 4 < (3$ < 0.07. ICaselhll2002l ). Thus, the gravitational 
instability scenario may be more common for the planet formation. 



4.3. Effects of Magnetic Field 



In this study, we ignored the magnetic field. The magnetic field plays an important role 
in both the early collapse and main accretion phases, and is related to the formation process of 
the star and circumstellar disk. For example, a large fraction of the cloud mass is blown away 
by the protostellar outflow that is driven by the Lorentz force from the circumstellar disk. 
Thus, the protostellar outflow is expected to be closely related to the star formation efficiency 
(IMatzner fc McKed |2000| ; iMachida et al.l 120091 ) . The p rotostellar outflow is also related to 



the determination of the circumstellar disk mass (e.g.. Ilnutsuka et al. Il2009i ). In addition, 
in the circumstellar dis k, the turbulence owing to the magneto-rotational instability (MRI; 
Balbus fc Hawleyll 19911 ) contributes to the angular momentum transfer that is related to the 



mass accretion rate onto the protostar. Moreover, the magnetic field supp ress fragmentation 



(Machida et al. 


2005bl: ] 


VIachi< 


a et al. 


2008; 


Hennebel 


e & Tevssier 


2008 


) and the formation 


of the disk ( 


Mellon & Li 


2008. 


2009: 


Duffin & Pudritz 


2009; 


Hennebelle & Ciardi 


2009 


)■ 



However, even when the molecular cloud is strongly magnetized, the mass of the first 
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core and protostar is, respectively, Mdisk ~ 0.01—0.1 M & and M ps ~ 10 3 M Q ( IMachida et al. 



20071 ). Thus, as in the case of an unmagnetized cloud, the circumstellar disk mass overwhelms 



the protostellar mass at the protostar formation epoch. This is because the mass and size of 
each object (first core and protosta r) are determined by the thermodynamics of the collapsing 
gas (IMasunaga fc Inutsuka I l2000l ). In the main accretion phase, however, the masses of 
the circumstellar disk and protostar in magnetized clouds may be different from those in 
unmagnetized clouds because the angular momentum is transferred by the magnetic effect 
such as magnetic braking, protostellar outflow and MRI. We will investigate the circumstellar 
disk formation in a magnetized cloud in detail in a subsequent paper. 



5. Summary 

To investigate the circumstellar disk formation and its properties, we calculated the 
evolution of the unmagnetized molecular clouds from the prestellar core phase until the end 
of the main accretion phase going through the protostar formation, using three-dimensional 
nested-grid simulation that covers both whole region of the initial molecular cloud and sub- 
AU structure. We constructed 10 models with a parameter of the initial cloud's rotational 
energy (/3o), in which the size (6.2 x 10 3 AU) and mass (1 M Q ) of the initial molecular cloud 
are fixed, and calculated their evolution until about 90% of the host cloud mass accretes 
onto either protostar or circumstellar disk (i.e., until the end of the main accretion phase). 
The following results are obtained: 

• The first core is a precursor of the circumstellar disk. In the collapsing cloud core, 
the first core forms at n c ~ 10 11 cm -3 with a size of ~ 10 AU before the protostar 
formation. At its formation, the first core exhibits a thick disk-like configuration with 
aspect ratio of H/R ~ 0.6 — 0.9. Then, after the protostar formation, the first core 
gradually transforms into a considerably thin disk (H/R < 0.1) that corresponds to 
the circumstellar disk. Thus, the minimum size of the circumstellar disk (or the size 
of the youngest circumstellar disk) is ~ 10 AU that corresponds to the size of the 
first core at its formation. The circumstellar disk is mainly supported by the thermal 
pressure gradient force in the early collapse and early main accretion phases, and then 
it is mainly supported by the centrifugal force and becomes very thin in the later main 
accretion phase because the gas with a large angular momentum accretes onto it. 

• The circumstellar disks are massive than the protostars throughout most of the main 
accretion phase. The first core has a mass of ~ 0.01 — 0.1 M Q , while the protostar has 
a mass of ~ 10 -3 M & at its formation. The mass of each object almost corresponds to 
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the Jeans mass at its formation. The first core evolves directly into the circumstellar 
disk. Thus, at the protostar formation epoch, the circumstellar disk is inevitably 
more massive than the protostar: the circumstellar disk is about 10 — 100 times more 
massive than the protostar. After the protostar formation (i.e., in the mass accretion 
phase), the gas accretes onto the protostar through the circumstellar disk. Since a 
part of the accreting gas accumulates in the disk, not only the protostellar mass but 
also the circumstellar mass increases. The mass growth rate of the protostar is larger 
than that of the circumstellar disk in the main accretion phase. Nevertheless, the 
circumstellar disk is more massive than the protostar throughout the main accretion 
phase, reflecting the large initial mass ratio between the protostar and circumstellar 
disk at their formation. 

• The planet and brown-dwarf mass objects tend to appear owing to the gravitational 
instability of the circumstellar disk in the main accretion phase. Since the circum- 
stellar disk is massive than the protostar throughout the main accretion phase, the 
gravitational instability inevitably occurs in such massive disk. In the cloud with ini- 
tially moderate rotational energy (10~ 5 < /3q < 5 x 10~ 3 ), fragmentation occurs in 
the circumstellar disk to form several clumps that have a mass of 20 — 70 Mj up with 
a separation of 7 — 100 AU from the protostar. On the other hand, in the cloud with 
initially larger rotational energy (5 x 10~ 3 < j3 < 3 x 10~ 2 ), fragmentation does not 
occur, but spiral patterns appear in the main accretion phase. Whether fragmentation 
or spiral arms appear depends on growth of the disk: fragmentation can occur when 
a sufficiently larger disk (rd; s k > A cr i) has sufficient mass (Q < 1). However, because 
the circumstellar disk not showing fragmentation in the main accretion phase has a 
mass comparable to the protostar, fragmentation may occur in later phase of the disk 
evolution, in which the disk cools and becomes more unstable again. 

• The mass accretion rate onto the protostar shows strong time variability in the main 
accretion phase. After low-mass (planet or brown-dwarf mass) objects appear in the 
circumstellar disk, the mass accretion rate onto the protostar shows a strong time vari- 
ability, because low-mass objects effectively transfer the angular momentum outwardly 
to promote gas accretion onto the protostar. On the other hand, there is little time 
variability of the accretion rate in a cloud with considerably lower rotation energy 
(flo = 10~ 5 ) in the initial state. In such a model, since neither clumps nor clear spiral 
patterns appear in the circumstellar disk, the angular momentum transport is not so 
effective in the present modeling without magnetic field. Thus, the variability of the 
accretion rate provides a useful signature for detecting the low-mass companions or 
exo-planets around very young protostars. 
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Our results imply that, in the star formation process, the planet and brown-dwarf 
mass objects frequently appear in the circumstellar disk, indicating that the gravitational 
instability scenario may be a major pass for the planet formation. To refine the formation 
scenario for planets, we need to study the disk evolution after the main accretion phase, in 
which both the radiative processes and magnetic effects should be taken into account. 

Numerical computations were carried out on NEC SX-9 at Center for Computational 
Astrophysics, CfCA, of National Astronomical Observatory of Japan, and NEC SX-8 at the 
Yukawa Institute Computer Facility. This work was supported by the Grants-in-Aid from 
MEXT (20540238, 21740136). 

REFERENCES 

Arquilla, R., & Goldsmith, P. F. 1986, ApJ, 303, 356 
Balbus, S. A., & Hawley, J. F. 1991, ApJ, 376, 214 

Bodenheimer P., Burkert A., Klein R. I., & Boss A. P., 2000, in Mannings V., Boss A. P., 
Russell S. S., eds, Protostars and Planets IV. Univ. Arizona Press, , p. 675 

Cameron, A. G. W. 1978, Moon and Planets, 18, 5 

Calvet, N., Hartmann, L., & Strom, S. E. 2000, Protostars and Planets IV, 377 
Caselli, P., Benson, P. J., Myers, P. C, & Tafalla, M. 2002, ApJ, 572, 238 
Duffin, D. F., & Pudritz, R. E. 2009, ApJ, 706, L46 

Durisen, R. H., Boss, A. P., Mayer, L., Nelson, A. F., Quinn, T., & Rice, W. K. M. 2007, 
Protostars and Planets V, 607 

Dutrey, A., Guilloteau, S., & Ho, P. 2007, Protostars and Planets V, 495 

Enoch, M. L., Corder, S., Dunham, M. M., & Duchene, G. 2009, arXiv:0910.27l5l 

Gammie, C. F. 2001, ApJ, 553, 174 

Goodman, A. A., Benson, P. J., Fuller, G. A., & Myers, P. C. 1993, ApJ, 406, 528 

Goodwin S. P., Kroupa P., Goodman A., & Burkert A., 2007, in Reipurth B., Jewitt D., 
Keil K., eds, Protostars and Planets V. Univ. Arizona Press, , p. 133 



- 19 - 

Hayashi, C, Nakazawa, K., & Nakagawa, Y. 1985, in Protostars and Planets II, ed. D. C. 
Black & M. S. Matthews (Tucson: Univ. Arizona Press), 110 

Hennebelle, P., & Teyssier, R. 2008, A&A, 477, 25 

Hennebelle, P., & Ciardi, A. 2009, A&A, 506, L29 

Hunter, C. 1977, ApJ, 218, 834 

Inutsuka, S., Machida, M. N., & Matsumoto, T. 2009, submitted flarXiv:0912.5439~j) 
Kalas, P., et al. 2008, Science, 322, 1345 

Kratter, K. M., Matzner, C. D., Krumholz, M. R., & Klein, R. I. 2009. larXiv:0907.3476l 
Larson, R. B., 1969, MNRAS, 145, 271. 

Machida, M. N., Matsumoto, T., Tomisaka, K., & Hanawa, T. 2005, MNRAS, 362, 369 

Machida, M. N., Matsumoto, T., Hanawa, T., & Tomisaka, K. 2005b, MNRAS, 362, 382 

Machida, M. N., Matsumoto, T., Hanawa, T., & Tomisaka, K. 2006, ApJ, 645, 1227 

Machida, M. N., Inutsuka, S., & Matsumoto, T. 2007, ApJ, 670, 1198 

Machida, M. N., Tomisaka, K., Matsumoto, T., & Inutsuka, S. 2008, ApJ, 677, 327 

Machida, M. N., Inutsuka, S., & Matsumoto, T. 2009, ApJ, 699, L157 

Mizuno, H., Nakazawa, K., & Hayashi, C. 1978, Progress of Theoretical Physics, 60, 699 

Matsumoto T., & Hanawa T., 2003, ApJ, 595, 913 

Marois, C, Macintosh, B., Barman, T., Zuckerman, B., Song, I., Patience, J., Lafreniere, 
D., & Doyon, R. 2008, Science, 322, 1348 

Masunaga, H., Miyama, S. M., & Inutsuka, S., 1998, ApJ, 495, 346 

Masunaga, H., & Inutsuka, S., 2000, ApJ, 531, 350 

Matzner, C. D., & McKee, C. F. 2000, ApJ, 545, 364 

Mellon, R. R., & Li, Z.-Y. 2008, ApJ, 681, 1356 

Mellon, R. R., & Li, Z.-Y. 2009, ApJ, 698, 922 



-20 - 

Meyer, M. R., Backman, D. E., Weinberger, A. J., & Wyatt, M. C. 2007, Protostars and 
Planets V, 573 

Natta, A., Grinin, V., & Mannings, V. 2000, Protostars and Planets IV, 559 

Neuhauser, R., Guenther, E. W., Wuchterl, G., Mugrauer, M., Bedalov, A., & Hauschildt, 
P. H. 2005, A&A, 435, L13 

Saigo, K., k Tomisaka, K. 2006, ApJ, 645, 381 

Shu, F. H. 1977, ApJ, 214, 488 

Perri, F., & Cameron, A. G. W. 1974, Icarus, 22, 416 
Toomre, A. 1964, ApJ, 139, 1217 
Vorobyov, E. I., & Basu, S. 2007, MNRAS, 381, 1009 
Vorobyov, E. I. 2009, ApJ, 704, 715 

Watson, A. M., Stapelfeldt, K. R., Wood, K., & Menard, F. 2007, Protostars and Planets 
V, 523 

Walch, S., Burkert, A., Whitworth, A., Naab, T., & Gritschneder, M. 2009, MNRAS, 400, 
13 



This preprint was prepared with the A AS IATgX macros v5.2. 



- 21 - 



Table 1: Model parameters and calculation results 



Model 


A> 
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[s- 1 




M ps [M ] 


M disk [M Q ] 


M rcd a [M ] 
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^frag 


[y] 


M frag [M ] 
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1 
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2.7 x 
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270 
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6 
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15 
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202 


1.5 x 


10 4 


0.037 
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0.74 
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8 
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14 


0.52 
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0.65 


117 


7.7 x 


10 3 


0.049 
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0.20 
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10 
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0.90 


0.08 


0.02 


0.10 
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b Residual mass [= (M ini - M ps - M disk )] 

a Disk-to-protostar mass ratio (— M^sk/ M ps ) 

c Fragmentation epoch after the protostar formation 
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Fig. 1. — Evolution of the circumstellar disk for model 6 (/3o = 10~ 3 ). The density (color 
scale and contours) and velocity (arrows) distributions on the equatorial (z = 0) plane are 
plotted. The spatial scale is different for each panel. The time elapsed after the protostar 
formation (t c ) and protostellar mass (M ps ) are shown in the upper part of each panel. 
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Fig. 2. — The density (color scale and contours) and velocity (arrows) distributions for the 
same model as in Fig. CD (model 6; 0q = 10~ 3 ) are plotted on the cross section of the y = 
plane at the same epochs as in Fig. [TJ 
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Fig. 3. — The radial distribution of (a) number density n c , (b) ratio of the radial to the 
azimuthal velocity R v (=v r /v^,), (c) radial velocity v r , and (d) azimuthal velocity (v^) at 
different epochs in model 6 (f3 = 1CT 3 ). In panel b, the velocity ratio R v inside the disk is 
not displayed to stress the disk surface. 
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Fig. 4. — The density distributions (color scale) for model 3 (/3q = 10~ 2 ) are plotted on the 
cross section of the z = (each upper panel) and y = (eac/i lower panel) planes. The grid 
scale of lower panels is eight times as large as that of upper panels. The time elapsed after 
the protostar formation (t c ) and protostellar mass (M ps ) are shown in the upper side of each 
panel. 
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Fig. 5. — Final states on the cross section of z = for models with Po — 10~ 2 (model 3), 
10~ 3 (model 6), 10 -4 (model 8), and 10 -5 (model 10). The density (color scale and contours) 
and velocity (arrows) distributions on the equatorial (z = 0) plane are plotted. The time 
elapsed after the protostar formation (t c ) and protostellar mass (M ps ) are shown in the upper 
side of each panel. The close-up view of the central region is plotted for models showing disk 
fragmentation. 
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Fig. 6. — The evolution of the disk mass as a function of time before (left) and after (right) 
the protostar formation for models with /3 = 10~ 2 (model 3), 10~ 3 (model 6), 10~ 4 (model 
8), and 1(T 5 (model 10). 
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Fig. 7. — The evolution of (a) residual mass, M res (= M ini — M disk — M ps ), (6) protostellar 
mass, M ps , (c) disk radius i?, and (d) aspect ratio of the disk H/R against time elapsed after 
the protostar formation for models with /3q = 10~ 2 (model 3), (model 6), 10~ 4 (model 
8), and 10~ 5 (model 10). 
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Fig. 8. — The evolution of the ratio of disk mass to protostellar mass against the protostellar 
mass for models with /3 = 10~ 2 (model 3), 10~ 3 (model 6), 10 -4 (model 8), and 10 -5 (model 
10). 
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Fig. 9. — The mass accretion rate (left axis) and disk mass (right axis) are plotted for 
models 3, 6, 8 and 10. Note highly time-dependent M ps driven by gravitational torque due 
to non-axisymmetric structure in the disk. 
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Fig. 10. — Averaged Toomre's Q (left panel) and the ratio of the critical Jeans length to 
disk radius (right panel) against the protostellar mass for models 3, 6, 8 and 10. 
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Fig. 11. — Final states on the cross section of z = (each top panel) and y = (each bottom 
panel) planes for models with /3q = 3 x 1CT 3 (model 4), 5 x 10~ 3 (model 3) and 3 x 1CT 2 
(model 1). The density distribution (color scale) is plotted in each panel. The time elapsed 
after the protostar formation (t c ) and protostellar mass (M ps ) are shown in the upper part 
of each panel. 



